Irreversibility, steady state, and non-equilibrium physics in relativistic heavy ion 

collisions 
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Heavy ion collisions at ultrarelativistic energies offer the opportunity to study the irreversibility 
of multiparticle processes. Together with the many-body decays of resonances, the multiparticle 
processes cause the system to evolve according to Prigogine's steady states rather than towards 
statistical equilibrium. These results are general and can be easily checked by any microscopic 
string-, transport-, or cascade model for heavy ion collisions. The absence of pure equilibrium states 
sheds light on the difficulties of thermal models in describing the yields and spectra of hadrons, 
especially mesons, in heavy ion collisions at bombarding energies above 10 GeV/nucleon. 
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The hypothesis that local equilibrium (LE) is attained 
by the system of two heavy ions colliding at ultra- 
relativistic energies is a basic assumption of macroscopic 
thermal- and hydrodynamical models of heavy ion colli- 
sions. The idea was pushed by Fermi [jjj and Landau [^]|| 
almost 50 years ago for hadron-hadron collisions. Despite 
the long history of theoretical and experimental attempts 
there is no rigorous proof of LE yet. The present paper 
shows that the irreversibility of multiparticle processes, 
proceeding e.g. via string decays, causes these systems to 
evolve according to Prigogine's Q steady state solution, 
rather than towards statistical equilibrium. 

Using Bogolyubov's hierarchy of relaxation times in 
non-equilibrium statistical mechanics [|| one usually con- 
siders the following scheme: Suppose that in the initial 
stage the system is far from equilibrium. To describe it 
one has to introduce a set of various many-particle dis- 
tribution functions rapidly varying in time. Then, due 
to interactions between the particles, correlations of the 
distribution functions occur within very short time in- 
tervals which are typically on the order of the collision 
time. 

This is the kinetic stage - all many-particle distribu- 
tion functions may be derived from the single one-particle 
distribution function. For times significantly larger than 
the collision time the number of parameters character- 
izing the system is reduced further to very few average 
values, namely the number of particles, their energy and 
velocity, i.e. to the moments of the distribution func- 
tion. At this stage the system behavior is governed by 
hydrodynamics. 

Unlike in non-relativistic mechanics, in relativistic 
heavy ion collisions the relaxation picture is more com- 
plex because of multiparticle processes. Here the num- 
ber of particles and their composition are not conserved. 
Newly produced particles are not thermalized (even if 
they appear to be, see and this circumstance causes 
a delay in achieving equilibration. The equilibration time 
may appear too long as compared to the typical lifetime 



of the expanding system. Due to the lack of a rigorous 
first-principles theory of nuclear reactions at relativistic 
energies, the approach to LE is investigated mainly by 
virtue of dynamical calculations provided by microscopic 
semiclassical Monte Carlo models PHl5|] which have been 
intensively studied during the last 15 years. 

The analysis of the space time evolution picture ob- 
tained in these models reveals that the whole system of 
colliding nuclei never attains a global equilibrium state 
after the initial non-equilibrium stage. Still, there is, in 
principle, a possibility of the occurrence of local equilib- 
rium (e.g. in the central cell), because the approach to 
LE does not depend on the assumptions of the presence 
of a heat reservoir, of Gibbs ensembles, etc. 

Our study has been inspired by the finding that quasi- 
stable states are present in partonic and hadronic mat- 
ter, as observed independently in dynamical simulations 
[p|, ^6|Jl7| . On the partonic level an analysis of the ther- 
malization of partons has been performed by the late 
Klaus Kinder-Geiger |8| . Equilibration of hadronic mat- 
ter has been studied, e.g., in the Quantum Molecular Dy- 
namics models ||{l(],[l3| |l5],[l8|. These simulations have 
shown that at high energies neither the global system nor 
its central part seem to reach the state of chemical equi- 
librium (in the sense of statistical mechanics) HQ . This 
observed feature is not solely restricted to microscopic 
models. To describe, for instance, the experimental data 
on yields of strange particles in heavy ion collisions at 
200 AGeV pi or hadron multiplicities at 158 AGeV pfj 



the standard statistical model of the ideal hadron gas 
has been modified to invoke the hypothesis of chemical 
non-equilibrium pT| , f22fl as well. 

Does this simply imply that the hadronization time is 
shorter than the equilibration time? - Not necessarily! 
In the present paper we show that dissipative processes, 
such as multiparticle production via strings and many- 
body (N > 3) decays of resonances, dominating at high 
energies, can lead to the creation of a stationary state 
called steady-state. This steady-state does not coincide 
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with a pure "conventional" equilibrium state, as assumed 
in the statistical models. 

Consider first the necessary and sufficient criteria of 
LE in the central zone of nuclear reactions, which is usu- 
ally analyzed in microscopic calculations: 
Necessary conditions: (i) absence of significant flow 
effect in the central cell; isotropy of the velocity distri- 
butions, and (ii) isotropy of the diagonal components of 
the pressure tensor, 
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Here V is the volume of the cell and rrii,Pi are the mass 
of the i-th hadron and its momentum, respectively. 
Sufficient conditions: (iii) thermal equilibration 
which manifests itself in the time independence of the 
hadronic spectra after a certain period, and (iv) chemi- 
cal equilibration, i.e. the time independence of different 
hadronic yields. 

The necessary conditions look quite simple and evi- 
dent: Local equilibrium may not be reached in symmet- 
ric nuclear collisions earlier than for the time t pass = 
2R/( , y cm v cm ), during which noninteracting Lorentz con- 
tracted nuclei of radius R/jcm, which stream freely with 
the velocity v cm , would have passed through each other. 
Apparently, early in the collision this is the origin of a 
substantial initial longitudinal collective flow of hadrons 
in the cell, which distorts the equilibration picture at 
the very beginning of the reaction. After t pass this non- 
equilibrium component rapidly drops |14|. In |l5[ ] it has 
been reported that a stage of kinetic equilibrium is at- 
tained in heavy ion collisions in a central cell of volume 
V = 5x5x5 — 125 fm 3 at about t = 10 fm/c, irre- 
spective of the energy of the colliding nuclei from 10.7 
AGeV (AGS) to 160 AGeV (SPS). Isotropy of both the 
pressure and the velocity distributions of hadrons charac- 
terizes, without the sufficient conditions (iii) and (iv), 
however, a pre-equilibrium stage of the reaction rather 
than an equilibrium one! In a fully equilibrated system 
conditions (iii)-(iv) must be satisfied as well. 

This is the crucial point in our discussion: The statisti- 
cal thermodynamics of many-particle systems |2^j2^| de- 
termines the thermal equilibrium as the state with max- 
imum entropy. Once thermal equilibrium is attained, 
the velocity distributions of different particles must be 
isotropic. If the total number of particles is conserved, 
kinetic equilibrium is equivalent to thermal equilibrium 
p4j . But: this equivalence is broken, both in chemical 
reactions and in high energy physics. 

Indeed, if the mixture of reacting substances is in the 
"true" equilibrium, then the rates of each chemical reac- 
tion must be the same for the direct and inverse processes 
[pj,^3|. However, in a cyclic process, in which the concen- 
trations of the reacting substances are time independent, 
but the partial reaction rates uij = uij" — uj™ are non- 



zero, the system is in a stationary state, which may be 
far from the equilibrium [EJ. 

Consider now an ideal thermostat which contains a few 
thousand protons with an energy E 3> m p c 2 , where m p 
is the mass of proton and c is the light velocity. For the 
sake of simplicity we exclude the (slow) weak processes 
from this scenario, focusing on strong interactions only. 
Then, even if the initial momentum distribution of the 
protons is Maxwellian, thermal equilibrium (in the sense 
of a state of maximum entropy) is not reached yet. Many 
new particles, mostly pions, will be produced as a result 
of initial proton-proton and, later, proton-pion, etc. col- 
lisions. When the system will finally reach equilibrium, 
it will consist of a large number of pions (and heavier 
mesons) with an admixture of baryons (and antibaryons) 
whose net number is conserved. The final temperature 
must, of course, be much lower than the initial one - 
kinetic energy has been transformed into mass (of pro- 
duced particles). But: will the particle abundances be 
the same as those given by the statistical mechanics of an 
ideal hadron gas? In other words, will the final state be 
the state of thermal and chemical equilibrium, in which 
any direct and inverse hadronic processes will be taking 
place on average at the same rate? 

This problem is closely related to the principle of de- 
tailed balance and to the irreversibility of multiparticle 
processes. To avoid ambiguities, we would like to stress 
that the definition of detailed balance in quantum me- 
chanics (DB^ M ) does not coincide with the definition of 
detailed balance in statistical physics (DB SP ) and chem- 
istry. Detailed balance in the sense of quantum mechan- 
ical invariance under time reversal implies that the tran- 
sition amplitudes of the direct and the inverse processes 
must be of the same magnitude, 
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In statistical physics and chemistry, the principle of de- 
tailed balance requires that (in thermostatic equilibrium 
of a system) every separate reaction between its compo- 
nents is in itself in equilibrium, i.e. the rates of the direct 
and inverse processes are the same. To clarify the differ- 
ence between DB^ M and DB SP , consider the process of 
multiparticle production, e.g. in string excitation: 
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According to Fermi's Golden Rule, the probability of n 
particle production reads 
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where p is the total four-momentum, pi is the four- 
momentum of i-th particle, and \M\ is the amplitude of 
the process. The last factor is the space factor, which is 
fully determined by the kinematics of the reaction [E5|. 
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FIG. 1. Irreversible process: recombination of sev- 
eral hadrons into two baryons. Schematic diagram of (a) 
many-particle collision of hadrons; (b) rearrangement of their 
quarks and subsequent annihilation of gq-pairs; (c) produc- 
tion of two quark-diquark strings, which shrink to (d) two 
baryons of high energy. 

Although \M a+b ^ n \ = \M a+b ^„\ and, therefore, DB QM 
is satisfied, the rates of the direct and the inverse pro- 
cesses, Ra+b^n and R a +b^n, are different, due to differ- 
ent space factors. This means that DB SP is not fulfilled. 
Note that the principle of detailed balance in particle 
physics has been verified for the reactions 

a + b < — ► c + d , (5) 
a + b < — ► c , (6) 

where a,b,c,d denote hadrons and their resonances. 
These processes are time reversible, because the space 
factors (or the densities of states) of the initial and final 
states are essentially the same. 

The space factors are rapidly varying functions of n. 
Therefore, the matrix elements \M\ may be replaced by 
average values, and a situation typical for statistical me- 
chanics is obtained: the probability of a state is pro- 
portional to the volume of the accessible phase space. 
In other words, multiparticle processes are irreversible 
in time, because they increase the local entropy of the 
system. Consequently, the processes of recombination of 
many hadrons into one string, and two strings colliding 
to form a couple of ground state hadrons (Fig. [l]) of high 
energy are strongly suppressed, because they violate (lo- 
cally) the Boltzmann if-theorem. 

The irreversibility of multiparticle processes (e.g., 
strings which provide a steady source of new particles) 
first drives and then keeps the hadronic system out of 
the total chemical equilibrium, i.e. out of the full de- 



tailed balance in the sense of statistical mechanics. On 
the other hand, the DB SP principle is the basic assump- 
tion of the statistical model (SM) of the ideal hadron gas 
p(i| p9[ and variations like the statistical bootstrap model 
(SBM) ||. Therefore, simply extracting the energy den- 
sity £, the baryon density ps and the strangeness density 
ps of the system at a given time and inserting these val- 
ues as an input into the statistical model will be giving 
misleading results until all multiparticle processes in the 
system will have ceased. 

Still, the conditions (iii)-(iv) may be fulfilled, even if 
the full detailed balance is not reached yet. Such states, 
which may be stable or not, but are out of local equi- 
librium, have been dubbed steady states of the system 
[Q. To decide whether or not a steady state is attained 
in a microscopic model of heavy-ion collision, the sys- 
tem must be compared with the quasi-equilibrated (in 
the sense of the criteria (i)-(iv)) infinite matter, as sim- 
ulated within the same microscopic model. In |fl5| it was 
shown that the yields and energy spectra of hadrons in a 
central cell are - after t = lOfm/c - very close to those 
values calculated for infinite hadron matter ]l6| with the 
same e, Pb and ps- This is a strong indication on the 
occurrence of a steady state. 

In conclusion, we have discussed the relaxation of 
hadronic matter produced in the central zone of heavy 
ion collisions in the energy range spanning from AGS to 
RHIC. Apparently, dissipative TV-body (N ^> 1) decays 
of strings and resonances, i.e. multiparticle processes, are 
irreversible in time: the probability of N particles 

1) to collide simultaneously in a small volume and 

2) to transform into a final state, which consists only of 
two particles of higher energies, 

drops extremely rapidly with rising N. 

Therefore, these processes drive the system towards a 
steady state. Due to the broken symmetry between the 
rates of direct and inverse processes, this steady state 
does not coincide with a pure equilibrium state. 

The conditions (iii)-(iv), often applied for the anal- 
ysis of local equilibrium, are generally weaker than the 
requirement of full local equilibrium usually imposed in 
the macroscopic models. At low energy densities, when 
multiparticle processes are rare, the steady state coin- 
cides practically with the equilibrium one. At higher en- 
ergy densities, the difference between the states becomes 
more and more significant. 

One characteristic feature of the steady state would be 
a strong enhancement of pions, accompanied by a sup- 
pression of (many-body decaying) resonances. This is 
due to the absence of an effective feeding mechanism. 
This feature of the steady state can explain the fact why 
conventional thermal models considerably underestimate 
yields of pions at energies of E > 10 AGeV. These results 
are typical for a large family of microscopic (cascade-, 
transport-, string-) models, which describe hadronic and 
nuclear interactions without invoking the hypothesis of 
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quark-gluon plasma (QGP) creation. 

Non-equilibrium thermodynamics of irreversible pro- 
cesses j| finally comes to high energy physics, where 
the conservation of mass and particle number, conven- 
tional in statistical physics, is obviously violated. The 
number of possible reaction channels is three order of 
magnitude higher than in simple chemical reactions (see, 
however, the role of the equilibrium concept in biochem- 
ical/biophysical processes). Therefore, it is a hopeless 
task to solve the rate equations analytically. On the 
other hand, microscopic models for hadronic and nuclear 
collisions provide a very useful tool to probe these fun- 
damental features of nature at very small space and time 
scale. The non-equilibrium aspects of heavy ion collisions 
are interesting and require further investigations. 
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